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From stellar spectra, a variety of physical properties of stars can be derived. In particular, the chemical com¬ 
position of stellar atmospheres can be inferred from absorption line analyses. These provide key information 
on large scales, such as the formation of our Galaxy, down to the small-scale nucleosynthesis processes that 
take place in stars and supernovae. By extending the observed wavelength range toward bluer wavelengths, 
we optimize such studies to also include critical absorption lines in metal-poor stars, and allow for studies of 
heavy elements {Z > 38) whose formation processes remain poorly constrained. In this context, spectrographs 
optimized for observing blue wavelength ranges are essential, since many absorption lines at redder wave¬ 
lengths are too weak to be detected in metal-poor stars. This means that some elements cannot be studied in 
the visual-redder regions, and important scientific tracers and science cases are lost. The present era of large 
public surveys will target millions of stars. It is therefore important that the next generation of spectrographs 
are designed such that they cover a wide wavelength range and can observe a large number of stars simulta¬ 
neously. Only then, we can gain the full information from stellar spectra, from both metal-poor to metal-rich 
ones, that will allow us to understand the aforementioned formation scenarios in greater detail. Here we 
describe the requirements driving the design of the forthcoming survey instrument 4MOST, a multi-object 
spectrograph commissioned for the ESO VISTA 4m-telescope. While 4MOST is also intended for studies of 
active galactic nuclei, baryonic acoustic oscillations, weak leasing, cosmological constants, supernovae and 
other transients, we focus here on high-density, wide-area survey of stars and the science that can be achieved 
with high-resolution stellar spectroscopy. Scientific and technical requirements that governed the design are 
described along with a thorough line blending analysis. For the high-resolution spectrograph, we hnd that a 
sampling of > 2.5 (pixels per resolving element), spectral resolution of 18 000 or higher, and a wavelength 
range covering 393-436 nm, is the most well-balanced solution for the instrument. A spectrograph with these 
characteristics will enable accurate abundance analysis (±0.1 dex) in the blue and allow us to confront the 
outlined scientihc questions. 
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1 Introduction 

Spectra obtained in the near-ultraviolet to blue wave¬ 
length range (300 ^ A ^ 450 nm) can carry a 
wealth of information. In stars, a blue colour can in¬ 
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dicate that the star is hot, metal-poor, or a com¬ 
bination of both. Most metals {Z > 3) have their 
strongest absorption lines in this wavelength region. 
The heavy elements {Z > 38) have their strongest, 
and sometimes their only absorption lines accessible 
from ground-based observatories, blue-w ards of 430 nm 
( Cowan et al. 2002 : ISneden et aL 20031) . This region is 
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typically very crowded causing many severe line blends. 
Therefore, this region is often avoided in spectral anal¬ 
yses, as they require detailed and accurate line lists to 
allow for a precise chemical abundance analysis from 
the atomic and molecular features. Furthermore, high- 
resolution spectra are also needed to better resolve 
and separate the individual absorption features of the 
blending components. 

Currently, there are only few high-resolution (i.e., 
R = A/AA ^ 20 000), blue-to-near-UV sensitive 
spectrographs on large ground-based telescopes (> 
4 m), namely the High-Resolution Echelle Spectro¬ 
graph (HIRES) at the 10-m Keck telescope, the High- 
Dispersion Spectrograph (HDS) at the Subaru 8.2-m 
telescope, and the Ulraviolet-to-Visual Echelle Spec¬ 
trograph (UVES) at the 8.2-m Very Large Telescope 
(VLT), all of which operate down to the atmospheric 
cut-off at ^ 300 nm. The MMT Advanced Echelle 
Spectrograph (MAESTRO) at the MMT0 covers the 
UV spectral range down to ^ 315 nm, and on the 
same telescope, the Multi-object Hectochelle down to 
380 nm. Finally, there are the Magellan Inamori Ky¬ 
ocera Echelle (MIKE) instrument at the 6.5-m Magel¬ 
lan telescope with a wavelength limit of 335 nm, and 
the High Resolution Spectrograph (HRS) at the 11-m 
South African Large Telescope (SALT) for which the 
bluemost observable wavelength is ^ 350 nm. 

However, in order to answer some of the important 
open questions in astrophysics related to the origin 
and formation of the Galaxy as well as the creation 
of heavy elements (see Sect. El), it is imperative to de¬ 
sign more blue-sensitive, high-resolution, multi-object 
spectrographs. In particular, in old, metal-poor stars as 
typically found in the Galactic halo, we need to study 
the blue wavelength ranges to detect some of the heavy 
elements, as their redder lines are often too weak to be 
detectable. To address these questions, we carry out de¬ 
tailed studies on when absorption lines are detectable 
and useful in the blue spectral range for an abundance 
analysis. 

The 4-metre Multi Object Spectroscopic Telescope 
(4MOST) will be a new fibre-fed instrument meant to 
simultaneously provide high- and low-resolution spec¬ 
tra. It is planned to be mounted at the VISTA telescope 
(Visible and Infrared Survey Telescope for Astronomy) 
at ESO’s Paranal Observatory. Gurren tly, 4MOST is 
in it s Prel iminary Design phase (see Ide Jong et al.l 
2 OI 4 L 2 OI 2 L for a detailed overview). It will carry out 


follow-up and extended observations for Gaia in ad- 
dition to the large ongoing Gaia-ESO Survey (GES; 
iGilmore et al. 2012 ). The high-resolution spectrograph 
design of 4MOST also includes a green and a red arm, 
which lead to a wide spectral coverage (from ^ 390 nm 
to ~ 676 nm, including signihcant gaps; Ruchti et al. 
2015, subm.), which allows for studies of substructure 


also known as Multiple Mirror Telescope 


in all major components of the Milky Way (MW), 
namely the disks, bulge, and halo. In addition, the in¬ 
strument will simultaneously obtain spectra at lower 
resolution {R > 5000), and will pursue both Galactic 
as well as extragalactic studies. Within the Galaxy, the 
large number of low- and high-resolution hbres will pro¬ 
vide spectra of up to 2400 stars spread over a field of 
at least 4 degrees^ field in one single pointing. Hence, 
4MOST envisions to target up to 30 million objects in 
its first survey of five year duration. 

Ruchti et al. (2015, subm.) describe how the red¬ 
der wavelength ranges (A > 450 nm) of the 4MOST 
high-resolution spectrograph were optimized to maxi¬ 
mize a versatile scientihc outcome. Here, we focus on 
the scientific requirements of the high-resolution blue 
arm that drives the design, how its wavelength cover¬ 
age is optimized for chemical tagging, and we describe 
blending issues that apply to every kind of spectrum 
densely populated by spectral lines. 

This paper is arranged as follows: Section[2]sketches 
a concept design and highlights potential construction 
issues. Section [3] outlines the scientific drivers. Sect. S] 
describes the scientific requirements such as wavelength 
coverage, resolution, sampling, and blending issues. 
Einally, the conclusion and outlook can be found in 
Sect. [3 

2 Constructing the blue arm of 
4MOST 



Fig. 1 Layout of the 4MOST HRS blue channel: 
common collimator optics, separation into three chan¬ 
nels by dichroics (here only the blue channel is shown), 
followed by the VPH grating, and finally the camera. 


The conceptual optical design and layout of the 
4MOST High-Resolution Spectrograph (HRS) is pre¬ 
sented in this section. We emphasize that the details of 
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Table 1 Conceptual design specifications of the 
4MOST HRS. 


Components/Parameters 

Values 

Wavelength range 

393-436 nm (blue) 
516-572 nm (green) 
610-676 nm (red) 

Spectral resolution 

> 18 000 (goal: 20 000) 

Fiber core diameter 

85 /rm 

Fiber separation 

170 /rm 

Fiber output focal ratio 

F3 

Spectrograph beam size 

250 mm 

Collimator 

Reversed off-axis 
Schmidt telescope 

Grating 

VPHG: 

3586 1/mm (blue) 

2836 1/mm (green) 
2343 1/mm (red) 

Camera focal ratio 

F1.8 

Detector 

6K X 6K, 15 ^m/pixel 

Sampling (spatial direction) 

3.4 pixels 

(dispersion direction) 

2.5 - 3.4 pixels 

Gap between fiber images 

3.4 pixels 

Number of fibres 

~ 840 


the design is conceptual, and the information provided 
is only intended to guide the reader. The HRS will 
have three channels covering the wavelength range from 
392.8 to 675.5 nm; we will focus here on the blue chan¬ 
nel covering the wavelength range from 393 to 436 nm. 

The main design parameters are given in Table [T] 
and the conceptual layout of the blue channel is shown 
in Fig.[T] Around 840 fibers can be injected in the HRS. 
After a common collimator optics, the light is split 
by dichroic filters into three separate channels to opti¬ 
mize the efficiency of the instrument. Volume-phased- 
holograpic (VPH) gratings are foreseen to disperse the 
light in each channel. The camera optics for all chan¬ 
nels is very similar, but not identical, and is also opti¬ 
mized for the wavelength range to be covered. As de¬ 
tector a CCD with 6kx 6k pixels is under evaluation. 
With this configuration, a minimum spectral resolving 
power of i? = 18 000 (goal: 20 000 - see also Sect. 14.1.111 
is reached over the entire spectral range covered. The 
theoretical optical performance of the system is excel¬ 
lent. The mean value of the width that receives 80% of 
the energy is 0.7 pixel over all wavelengths and fields. 
The poorest value is obtained at the CCD edge and 
corresponds to 1.5 pixels. Thus, minimal degradation 
of the nominal resolution/spectral purity will occur. 

In the blue wavelength range, the efficiency is a 
point that needs special attention. The total efficiency 
of the HRS blue channel, not including atmosphere and 
telescope, is given in Table [2j For a fiber-fed spectro¬ 
graph such as 4MOST, the major limitation comes from 
the fiber throughput. Therefore, special measures are 
taken in the layout of the instrument at the telescope, 


Table 2 Total efficiency of the blue channel not in¬ 
cluding effects from the atmosphere and telescope. 


Blue channel 

Component 

393 nm 

414 nm 

436 nm 

Fibers/Coupling 

0.74 

0.77 

0.79 

Collimator 

0.95 

0.95 

0.95 

Dichroic(s) 

0.98 

0.98 

0.98 

VPH grating 

0.49 

0.75 

0.48 

Camera 

0.92 

0.93 

0.94 

Detector 

0.88 

0.92 

0.91 

Total 

0.27 

0.46 

0.30 


such as to minimize the length of the fibers, which 
for 4MOST will be about 20 m. In Fig. [21 the internal 
transmission of a 4MOST fibre candidate for a length 
of 20 m is shown with the HRS wavelength ranges in¬ 
dicated. While the throughput is well above 90% for 
the green and red ranges, it drops to a mean value of 
about 85 % in the blue range. 



300 400 500 600 700 

Wavelength [nm] 


Fig. 2 Expected throughput for a fibre candidate to 
be used in 4MOST with a length of 20 m. The curve 
shown is only internal transmission, no coupling losses 
are included. The boxes indicate the wavelength ranges 
for blue the (A ^ 400nm), green (A > 500nm), and red 
(A > 600 nm) channel, respectively. 

A typical, indicative detector quantum efficiency is 
shown in Fig. [3] There are several very efficient coat¬ 
ings for the detector, giving a high efficiency over the 
full HRS wavelength range: a coating candidate is illus¬ 
trated in Fig. |3l which in the blue yields an efficiency 
close to 90 %. For the VPH grating, the theoretical ef¬ 
ficiency was calculated and a contingency factor of 0.9 
applied, yielding a conservative estimate. The through¬ 
puts of the common collimator, dichroics and camera 
are summarized in Table |2j A peak efficiency of 46 % 
is expected, with 27%/30% at the ends of the blue 
wavelength range. 
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Wovelength [nm] 

Fig. 3 Indicative detector efficiency for 4MOST HRS 
on the basis of a typical 6 k x 6 k chip and a candidate 
coating. The boxes indicate the wavelength ranges (as 
in Fig. [2|) for the blue, green, and red channel, respec¬ 
tively. 


3 Science drivers for the blue region 

The Milky Way (MW) provides a unique study case for 
the formation of barred spiral galaxies. Detailed infor¬ 
mation on the formation history can only be accurately 
extracted from stars in our Galaxy, constraining mod¬ 
els of spiral galaxy formation. The understanding of the 
involved physical processes can then be generalized to 
other galaxies that are too distant to be studied with 
the same accuracy. 

Of particular importance in such studies are the 
chemical abundances of individual stars. For exam¬ 
ple, the a-elements (such as O, Mg, Si, Ca, and Ti) 
provide information on whether or not the star was 
likely formed in situ or brought into the Galaxy af- 
ter a merger event with dwarf ga l axy-li k e objects (e.g., 
Kochi 2009: Nissen fc Schusteil I 2 OIOI : ISearle fc Zinn 


19781 : I Venn et al.l 120041) . The general low star forma¬ 


tion efficiencies in the low-mass dwarf galaxies and the 
ensuing paucity of the a-producing Type II Supernovae 
(SNe) at the time the Fe-producing S Ne la start to con¬ 
tribute leads to a l o w g/Fe-ratio (e - g., Aden et ^l201ll : 
iKoch et ^ 2008 : Shetrone et al.l I2OO3I) . in contras; 
to the rapid enrichment the halo experienced creat¬ 
ing higher [g/Fe] a,bund ances fe.g.. iNissen fc Schuster 
2OIOI : iRuchti et al. 2014 ). 


At low metallici ty (;$ —2 dex) the low-luminosity 
dwar f galaxies (e.g., iGilmore et al.ll2013l : iKoch fc Rich 
2 OI 4 I) share the same enhancement in a-elements as is 


kinematical analyses. Thus tracing, e.g., stellar streams 
from disrupted satellites in phase space, ancient disrup¬ 
tion events that are still coherent in velocity space can 
be efficiently recovered. In turn, these accurate kine¬ 
matic tracers, such as streams and stellar velocities, 
are a requirement to constrain the mass distribution of 
the Milky Way and the shape of its dark matt er halo 
re.g.. iKopospy et 20 lot Williams et al. 2011 ). 

Altogether, this will also enable a characterization 
of the progenitor systems such as their mass, lumi¬ 
nosity, star formation rate, and chemical enrichment 
histories. In particular old, metal-poor stars will pro- 
yide a better understanding of the early formation 
and chemical eyolution of the Galaxy. This outlines 
the importance of accurate dynamics and abundance 
studies, for whic h high-resolution spectrogr aphs are 
needed (see also Lindegren fc Feltzingl l2013l) . Hence, 
high-resolution stellar spectra yielding radial yelocities 
and elemental abundances will characterize a yariety 
of features from the nuclear formation processes that 
created the elements to shed light on large-scale galaxy 
formation scenarios. 

The stellar halo is the largest (yet least massiye) 
component in the MW. There is growing eyidence from 
kinematical and chemical abundance studies that it 
consists of two distinct subcomponents, an inner and 


2007: Nissen fc Schuster 


an ou ter halo (|Garollo et al 

2 OIOI). although these results are still being debated 


seen in Galactic halo stars (e.g., ICayrel et al. 2004 ). 
Howeyer, this is not the case at higher metallicities, 
where the stars in the dwarf galaxies tend to show 
lower a-abundances compared to the Galactic stars. 
Furthermore, understanding the assembly history of 
our Galaxy also benefits from mapping the dynam¬ 
ics of important Galactic substructures through stellar 


(|Sch6nrich et al.l|2011 ). These two components are sep¬ 
arated in their kinematics and elemental abundances, 
in particular iron abundances, from which a metallic- 
ity distribution function (MDF) can be calculatecd. 
Studies haye shown that the outer halo is charac¬ 
terised by metal-poor stars (< [Fe/H] —2.2) on 
retrograde orbits, indicating that these stars were ac¬ 
creted from other galaxies, while the inner halo stars 
(< [Fe/H] >~ —1.6) s how prograde orbit s, implying 
that they formed in situ ( Garollo et al.ll2007 ). A deeper, 
large survey of stars in both the inner and outer halos 
is needed to improve our current, incomplete under¬ 
standing of hierarchical galaxy formation, and metal- 
poor stars play a crucial role for understanding this. 
4MOST will consist of low- {R ^ 5000 with a goal 
of 7500) and high-resolution {R ~ 18 000; goal 20 000) 
spectrographs that are optimal for finding metal-poor 
stars with B > 17.5 mag, and acquiring high-resolution 
spectra of promising metal-poor candidates down to 
B = 16 mag. 

In general, the halo is thought to be the Galac¬ 
tic component containing the oldest stars, and they 
are metal-poor. Recently, also metal-poo r stars in the 
direc tion of the bulge have been found ( Howes et al.l 
2014I). and its old age has been demonstrated in 
Clarkson et alJ ( 2008[) . The spectra of these inetal- 


^ This MDF agrees with the one obtained from low-resolution 
metallicity indicators such as the Ca triplet or the Ca HK lines. 
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poor stars show very few detectable m etal lines 
(jChristlieb et alJl2004 iHansen et alJl2014br) . and some 
of the strongest lines, still detectable in metal-poor 
stars, are found in the blue spectral region. In order to 
extract the maximum possible information from these 
metal-poor halo stars, a blue wavelength range is es¬ 
sential. 

This is particularly true for the absorption lines of 
the heavy elem ents, as most of these lines are found be - 
low ~ 430nm ( Cowan et al. 2002 : Sneden et al. 2003). 


They provide important tracers of elements created 
by the neutron-capture processes. There are two main 
channels for creating neutron-capture elements; the 
rapid and slow neutron-capture process. The exact 
physics and sites of these processes are poorly con¬ 
strained, and this emphasises the need for abundances 
of heavy elements from blue high-resolution spec- 


Hansen et al. 

2014allbl. 2012: 

Sneden et al. 

l2008l). 


Additionally, the blue range contains numerous el¬ 
ements useful for chemical abundance analyses of stars 
with [Fe/H] < — 1, but this range can also be used for 
chemical analyses of more metal-rich stars. The blue 
spectral range hosts not only transitions of neutron- 
capture elements, but also lines from elements be¬ 
longing to each of the following groups: a, odd-Z, 
(in)complete Si-burning, and Fe-peak elements, mak¬ 
ing this a crucial wavelength region for chemical tag¬ 
ging. Furthermore, the blue range A < 450 nm con¬ 
tains the A^A — X^Il G-band of CH, which can be 
used to detect carbon-enhanced metal-poor (CEMP) 
startH, and to determine their carbon abundances. 
Combined with the s-process elements (e.g., Sr, Y, 
Zr, Ba, La) and r-process elements (e.g., Eu, Nd), 
the blue spectra thus serve to classify CEMP stars 
into —s,—r,—rs, and —no subg roups through their 


heavy element a.bund ances (e.g., iBisterzo et al. 2012 : 


(20-40% or more: Beers & ChristliebI 20051: 

Lee et al. 

20131 iLucatello et al. 200fil Placco et al. 

20141 of 

the extremely and ultra metal-poor stars are en¬ 
hanced in carbon, and almost all the ultra Ee-poor 
stars are enhanced in carbon (^Bonifacio et al. 2015: 

Christlieb et alJ 2004: Frebel et al. 12007 

: Hansen et al. 

2014(1 Keller et al. 20141 Norris et al.l 

20071. Under- 


standing how these stars form is key to understanding 
how the early halo was enhanced in metals, and this in 
turn will place constraints on the physics of the forma¬ 
tion of the first few generations of stars in the Universe. 

All of the above underlines the need for spectro¬ 
graphs sensitive at blue wavelengths (A 450nm). 


^ M etal-poor stars with [C/Fe]> 0.7 according to llAoki et al.l 
l2007li . 


and 4MOST would only be the eighth spectrometrc0 
sensitive in the blue. 


4MOST also envisions a low-resolution mode {R = 
5000-7500) to cover a broad range of Galactic science 
goals. Many of these differ from the high-resolution sci¬ 
ence cases, so we will only briefly comment on the sci¬ 
ence facilitated by the blue region of the low-resolution 
setting. 


Even though the blue region can be used for chemi¬ 
cal tagging of more metal-rich disk stars, this region is 
suboptimal owing to the extreme line blending in the 
blue region in metal-rich stars (see Sect. 14.3.11 - 14.3.21 
for details on line blending). Therefore, chemical abun¬ 
dance analyses of metal-rich Population I stars are gen¬ 
erally conducted at redder wavelengths. Moreover, ra¬ 
dial velocities can also be determined from blue spectra 
of metal-rich disk stars, but again due to the large num¬ 
ber of heavily blended and unresolvable lines, leading 
to more uncertain measurements, the redder regions of 
the low-resolution mode of 4MOST, such as the region 
around the near-infrared Ca triplet lines, are preferred 
for this purpose. 


An important driver for the blue, low-resolution 
spectrograph is the detection and first characterization 
of hyper metal-poor star^ in the MW halo. Here, the 
blue region of choice (goal) comprises 390-543 nm, to 
include the Ca H and K lines (at A = 397 nm and 
393 nm, respectively) as the strongest lines that re¬ 
main detectable at these low metallicities. The Ba II 
455.4 nm line is an important tracer of neutron-capture 
processes, in particular of the earliest nucleosynthesis 
in the Universe, as it remains measurable down to ex¬ 
tremely low metallicity in most stars at the low reso¬ 
lution planned for 4MOST. This line will help us sub¬ 
classify the CEMP stars and, in turn, place constraints 
on their formation scenario. The red wavelength cut of 
this region is set by the Mg b triplet and the molecular 
CN-bands. Depending on the effective temperature and 
carbon-enhancement, also the G-band of CH (431 nm) 
remains detectable; it can be used for inferring C abun¬ 
dances to a precision of up to ±0.2 dex. Note that the 
Ca H and K lines can also serve as a measure of stellar 
activity in young stars. 


Only UVES is accessible to the entire astronomy commu- 
nity, while all other high-resolution spectrographs listed above 
are restricted-user access instruments. This would make 4MOST 
the second instrument with a blue arm that will be globally 
available. 

^ Following the definitions in iBeers & Christli^ 1120051) these 
stars have [Fe/H] < —4. 
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4 Scientific requirements to the blue 
instrument 

4.1 Sampling and resolution 

The spectral resolving power, R = A/AA, and sampling 
s (pixels per resolving element) are amongst the most 
important quantities that drive the design of a spectro¬ 
graph, and the science that can be achieved. In order to 
determine optimal values, taking into account technical 
as well as financial boundary conditions, 
an ap proach similar to that outlined in 
( 2013ll . In the following we will only give a brief sum¬ 
mary of our procedure, while our detailed tests will be 
presented in a forthcoming paper. 


4.1.1 Synthetic lines and parameter space 

First, we synthesized a Gaussian line profile, cen¬ 
tered at Ao, of given equivalent width (EW) and Full 
Width at Half Maximum (FWHM) at initially very 
high sampling over a continuum level of unity. The 
FWHM is chosen to be the resolution element (RE) 
of the spectrograph, so that the resolving power of 
the simulated spectrum can then be expressed as R = 
Aq / RE = Aq / FWHM [nm]. Next, the sampling of the 
simulated spectrum is degraded to a varying number of 
pixels per RE, and Poisson noise is injected to mimic 
different signal-to-noise ratios {S/N). The goal is to 
find a minimum set of requirements, that will still al¬ 
low us to derive abundances with an uncertainty of 
±0.1 dex or better. This corresponds to a relative EW 
uncertainty of 0.22-0. 2^. 

Eor testing purposes, we ran 10 000 Monte Carlo re¬ 
alizations with a fixed S/N of the 256 different combina¬ 
tions of each of the parameters in steps of: EW = (1.25, 
2.5, 5, 10) pm; i?= (15 000 - 23000); s = (l - 4) pix- 
els/RE; and S/N = (50, 60, 75, and 100 per pixel, or 
154, 185, 231, 308 per This test was performed to 
explore the parameter space, and following find a set of 
useful parameters to explore in more detail. In each of 
the Monte Carlo trials, the position of the resampled 
pixels was rigidly shifted by a random (linear) amount 
of up to ±0.5 times a resampled pixel with respect to 
the default value. This is to ensure that the position of 
the line center is different each time with respect to the 
positions of the resampled pixel, so as to avoid biases 
due to constant line-pixel offsets. 


® for weak lines and less for strong or partly saturated lines, 
we used S/Npi^^i oc S/N a y/ and we calculated the S/N 

per A for this specific example with a fixed value of s = 2.1, it = 
19 000, and A = 4200A . 


4.1.2 Line measurement, optimal resolution 
and sampling 

The resampled and noise injected lines were evaluated 
by fitting a Gaussian line profile within ±3xRE around 
the nominal line center. The continuum level, in turn, 
was estimated in two windows of 5xREs width and 
located at ±4xRE from the line center. 

The first tests showed that a sampling of 1 pixel/RE 
is simply too low to retrieve abundances within the 
required accuracy, and we find that the sampling needs 
to be between the lower limit of 2 (Nyquist sampling) 
and 3 pixel/RE. 

After exploring the parameter space, a more de¬ 
tailed test was carried out on the constrained set of 
promising parameters using finer steps compared to 
the aforementioned first test. Eigure[3]shows our results 
(using the restricted sampling of s = 2 and 3 pixel/RE) 
in terms of the relative EW measurement uncertainty, 
crEW/EW, as a function of the tested resolution. A line 
of 25 mA was chosen to mimic a weak line such as the 
670.6 nm Li line. 



15000 16000 17000 18000 19000 20000 21000 22000 23000 

Resolution 


Fig. 4 Relative EW uncertainty vs. resolution for 
different rates of S/N ratio, and a sampling of 2 and 
3 pixel/RE for the filled and open circles, respectively. 
An indicative relative EW uncertainty is plotted as a 
dashed, grey line. 

Imposing our requirement that lines need to be 
measured to a precision of ±0.1 dex elemental abun¬ 
dance to perform meaningful science for Galactic ar¬ 
chaeology, and that the S/N should be obtainable 
within a reasonable exposure time (e.g., S/N=170 per 
A or 55 per pixel), a trade off must be made. In con¬ 
junction with the preliminary instrument design, we 
chose the value of 2.5 pixel/RE as a compromisqj. This 
leads to our final specifications, namely: 

® We set a goal of 3 pixel/RE instead of the Nyqist-optimal of 
2 pixel/RE to maintain resilience against CCD cosmetics defects 
that may fall over a measured line. Dispersion variations along 
the spectral range can make sampling locally worse, so here we 
assume 2.5 pixel/RE as a compromise. 


we followed 
Caffau et al. 
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— R > 18 000 everywhere in the covered wavelength 
range; 

— a sampling of at least 2.5pixel/RE (goal: 
3.0pixels/RE). 


The adopted (goal) abundance uncertainty of 
0.1 dex is a value that, provided we have accurate 
atomic physics input, is well documented. It arises as a 
typical uncertainty owing to the influence from the stel¬ 
lar parameter uncertainty. This has been t ested in great 


detai l for 18 elements (see Table 6 and 7 in ICavrel et al 


20041 1. In Sect. 14.3)^ we will investigate how this uncer¬ 


tainty affects the usefulness of blends. 


4.2 Selecting and optimizing the blue 
wavelength range 


In order to fulfil the science requirements for the 
4MOST high-resolution survey focussing on the halo, 
elements from each major group (a-, odd-Z, iron-peak, 
and neutron-capture elements) should be included in 
the covered wavelength range. Such elements can be 
efficiently measured from absorption lines in the blue 
range of metal-poor halo stars, whereas in more metal- 
rich stars from the disks and bulge, a redder wavelength 
region is preferred owing to fewer line blends (Ruchti et 
al. 2015, subm.). As part of the 4MOST requirements, 
we want to be able to derive accurate abundances 
(O.Idex uncertainty) from several elements belonging 
to each group. Here we test the 390-460 nm (goal) re¬ 
gion and require that each absorption line reaches 0.9 
in residual flux, to be detectable even in noisy spectra. 


Synthetic spectra we re created wi th the MOOG 
synthetic spectrum code (|Snedenlll973l version 2010), 
in conju nction with the MARC S stellar atmosphere 
models ( Gustafsson et al. 20081) matching the stel¬ 
lar parameters listed below. A line list downl oaded 
from the VALD database i Kuoka et al.l l2000l) was 
used for this purpose. Several transitions in the line 
list were checked a nd updated with th e atomic data 
given in NIS T, and Cowan et alJ (12002 ). Sneden et al 


( 2003 ). and McWilliam et al. ()l995l) . The two first 
mentioned studies foc us on heavy elements, while 
iMcWilliam et al. ( I995h provides lists of both light and 
heavy elements and shows the detectability of numer¬ 
ous elements observed in various stars at a resolu¬ 
tion matching that of the 4MOST HRS. The goal is 
to obtain a set of useful, detectable absorption lines 
of: Mg, Al, Si, Ca, Mn, Fe, Co, Ni, Sr, Y, Zr, Ba, 
and Eu. To probe the detectability, we inspected syn¬ 
thetic spectra in which all abundances are scaled to 
Solar (i.e., [element/Fe] = 0.0). We target lines that 
are as clean as possible, reaching the aforementioned 
0.9 in residual flux, and we avoid saturated lines, to 
ensure that the lines are suited for chemical abun¬ 
dance analyses. This means that we need at least 
a S/N of about 30 (per pixel) to ensure Scr detec- 
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- Col 
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[Fe/H] 

Fig. 5 Detectable lines in the blue wavelength range 
in dwarfs (5500K/4.0/-I,-2,-3/2.0kms“^ - dashed, 
black line) and giants (4500K/2.0/-I,-2,-3/2.0kms“^ 
- solid, red line). Only a few Fe lines have been shown 
for simplicity. The yellow region indicates lines that 
could be lost, if the star has a strong G-band. 
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tions. In Fig. [SJ we show the detectability in a typi¬ 
cal dwarf and giant star as a function of metallicity, 
where the stellar parameters (Teff/log(?/[Fe/H]/Vt = 
5500/4.0/-1, -2, -3/2.0kms-\ and 4500/2.0/ - 1, 

—2, — 3/2.0kms“^, respectively) have been tested. 
T hese stellar paramet ers are consistent with the study 


of Caffau et al. ( 2013ll . 


Furthermore, also the radial velocity needs to be 
known to better than 1 — 2kms“^, a value that was 
set to be able to accurately trace kinematic substruc¬ 
tures such as streams from disrupted satellites in phase 
space. In this way, ancient disruption events that are 
still coherent in velocity space can be efficiently recov¬ 
ered. Finally, as we will see in Sect. 14.3.21 an accurate 
radial velocity is also important to recover abundances 
for blended lines. 


For finding the optimal wavelength range, a trade 
off between an as large as possible number of included 
lines on one hand, and fibre efficiency and the costs of 
the blue arm of the instrument on the other hand, has 
to be considered. Taking these constraints into account, 
we arrived at the wavelength region 393-436 nm (af¬ 
ter scanning the entire goal region: 390-460nm). This 
range ensures that at least one element of each ele¬ 
ment group is detectable, and can be measured in ei¬ 
ther dwarfs, giants, or both, at the required 4MOST 
resolving power of i? > 18 000. 

A few elements (Ni, and Ba) cannot be measured at 
this resolution in this rangcj^, and we need to use lines 
from the green and red wavelength ranges (Ruchti et 
al. 2015, subm.). However, in addition to the elements 
listed above, the blue range also offers a large number 
of other elements such as C (CH), Sc, Ti, Cr, La, Nd, 
Ce (illustrated in Fig. [5]). 

With a considerable number of both lighter and 
heavy elements well represented in the blue wavelength 
range, we can use the 4MOST spectra to accurately 
map the fossil gas in old stars. Our optimized range will 
also provide a considerable number of Fe and Ti I and 
II lines (minimum 7 and 5, respectively, covering exci¬ 
tation potentials of ^ 1 —3eV), needed for determining 
effective temperatures via excitation equilibrium, and 
surface gravities via ionisation equilibrium. Only the 
cleanest Ti and Fe lines are shown in Fig. [S] for the 
sake of clarity. This means that the blue region can, on 
its own, provide sufficient information needed for a full 
spectrum analysis. This is especially the case in metal- 
poor stars at [Fe/H] ^ —I, and the total number of 
clean lines will increase when including the green and 
red wavelength ranges, as envisioned for 4MOST. 

The HRS blue wavelength range was also set such 
that key elements like Al, Y, Sr, and CH (plus contin¬ 
uum redward of the band head), which are located at 
the edges of the wavelength range, would remain inside 


® That is, unless the star is enhanced in either of the elements 
or if the range is extended to 460 nm (goal region). 


Table 3 List of lines (elements, Z) tested for de¬ 
tectability in dwarfs/giants at [Fe/H]=—1, —2, —3. The 
‘worst’ losses are indicated with Columns three and 
four fall outside the covered wavelength range. 


Wavelength [nm] 

Z 

Wavelength [nm] 

Z 

394.4 

Al I 

439.8 

Y II 

395.0 

Y II 

442.9 

La II 

396.1 

Al I 

443.5 

Eu II 

397.1 

Eu II 

444.3 

Ti H 

398.8 

La II 

444.6 

Nd II 

399.5 

Co I 

444.9 

Dy H 

399.5 

La II 

445.4 

Ca I 

399.8 

Ti I 

445.9 

Ni 1* 

399.9 

Ce II 

446.1 

Fe I 

400.5 

Fe I 

446.2 

Nd II 

400.5 

V I 

446.8 

Ti II 

401.2 

Ti II 

448.2 

Fe I 

401.9 

Th II 

449.4 

Fe I 

402.1 

Nd II 

450.1 

Ti II 

402.3 

Nd II 

453.4 

Ti II 

402.8 

Ti II 

455.4 

Ba ID 

403.0 

Mn I 

456.3 

Ti II 

403.3 

Mn I 



403.4 

Mn I 



405.0 

Zr II 



405.7 

Pb II 



406.1 

Nd II 



407.7 

Sr II 



407.7 

Dy II 



408.6 

Th II 



408.6 

La II 



409.0 

Zr II 



410.9 

Nd II 



412.1 

Co I 



412.3 

La II 



412.9 

Ba II 



412.9 

Eu II 



413.7 

Ce II 



414.3 

Fe I 



416.7 

Mgl 



419.9 

Fe I 



420.2 

Fe I 



420.5 

Eu II 



420.8 

Zr II 



421.5 

Sr II 



422.6 

Ca I 



423.3 

Fe II 



424.6 

Sc II 



425.4 

Cr I 



427.4 

Cr I 



428.9 

Cr I 



429.0 

Ti II 



430.0 

Ti II 



430.2 

Ti I 



430.6 

Ti I 



430.8 

CH 



431.7 

Zr II 



431.8 

Ca I 



431.8 

Sm II 



432.2 

La II 



432.9 

Sm II 



433.3 

La II 
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the observed range even for high radial velocity stars 
which are common in the Milky Way halo. A wave¬ 
length piece of ^ 20 A is required in the beginning and 
end of the wavelength range before/after the first/last 
line of interest. This allows for loosing 5-10 A in the be¬ 
ginning of the spectrum to noise and even in stars with 
radial velocities of ^ 1000km s“^, we will still be able 
to get clean detections of A1 in their stellar spectra. 

The lines checked in the blue range can be found 
in Table [H Several of the weak lines may be detectable 
or useful if the star is enhanced in the element, or if 
the S/N ratio is larger than 30-50 per pixel. Not all 
the weak lines have been listed, as the line list contains 
more than 4000 lines, most of which are too weak to be 
detected in several of the tested stellar spectra at this 
resolution and S/N ratio. 

4.3 Blending issues 

The blue part of the spectrum is a very crowded re¬ 
gion, and as such, most absorption lines blend into 
each other. To investigate how strongly blended two 
lines can be, and still be useful for abundance de¬ 
termination, we start by investigating the idealized 
problem of htting two lines with Gaussian profiles, 
G(A) ~ a ■ exp A The goal is to obtain a 

relation between the line separation and the EW uncer¬ 
tainty resulting from the mutual blending. Following, 
we test this theory using synthetic spectra mimicking 
those from the future 4MOST HRS. 

4.3.1 A simple model of a two line blend 

In this and the following subsection, we consider two 
weak lines (a, b) with Gaussian line profiles, and we 
assume that there are no other blending lines interfer¬ 
ing with our tests. Furthermore, we assume that the 
stellar parameters, atomic data, and radial velocity are 
perfectly constrained, so that the line centres are at the 
rest wavelength. We want to find the minimum line sep¬ 
aration needed to fully recover the abundances of two 
blended lines. 

We consider two lines in a normalized spectrum 
characterized by two profile functions A > 0 and B > 0, 
scaled by amplitudes a and 6, respectively. We assume 
that the spectrum is sampled at m points, and that 
the noise in each pixel is independently normally dis¬ 
tributed with inverse S/N ratio n. We represent the 
lines in “depression format”, i.e. fitting the flux dif¬ 
ference to the continuum. The data model F of the 
spectrum is in this case 

F{X) = aA{X) + bB{X) (1) 

where A is the wavelength. We want to determine the 
amplitudes (a, b) scaling the line profiles of the ob¬ 
served data to be htted /,. Later, we specifically as¬ 
sume Gaussian profile functions {A, B) of given widths 


{wajWb) at wavelength positions (Aq, Xb). We determine 
a and b by maximum likelihood (hereafter ML) estima¬ 
tion. The problem is linear in a and b and we obtain 
as condition for a at the maximum of the (log) likeli¬ 
hood L 



( 2 ) 


with the analogous relation for b. 

We determine the uncertainties of a and b in the 
standard way from the inverse of the Fisher information 
matrix of the problem and obtain 


<> 2 - 




- B,A,Y 


(3) 


and 


(4) 


As expected, the standard deviations scale inversely 
proportional to the S/N. Strictly speaking, relations |3] 
and m give only lower (the Gramer-Rao) bounds of the 
variances of a and b. However, ML estimators are effi¬ 
cient and usually get close to the Gramer-Rao bound, 
so that we consider the lower bounds as actual uncer¬ 
tainties. 

We are interested in how the S/N ratio influences 
the accuracy to which the individual amplitudes a and 
b can be measured to ultimately obtain an estimate of 
the associated uncertainties in the derived elemental 
abundances. The problem of determining the EW of 
a blended line can be separated into determining the 
isolated clean line’s uncertainty stemming from noise, 
and finding a geometrical factor describing the blending 
configuration. When fitting a single line in isolation 
our simple model gives an uncertainty of the amplitude 
a' (assuming a profile function A) 


2 \ 2 
a„, > n 




(5) 


To have simple analytic expressions at hand, we now 
approximate the sums in the variances by integrals as¬ 
suming Gaussian line profiles of a given, unique width 
(iCa = Wh = w). For that, we have to assume a suf¬ 
ficiently fine sampling of the line profile, similarly to 
what was do ne in the derivation of Gayrel’s formula 
( GavrellflQSSh . We obtain for the uncertainties 0 -^ 
cr^, the relation 



( 6 ) 


Interestingly, equation (|6]) states that the detrimen¬ 
tal effect of blending is independent of the strength of 
the blending line as long as it remains on the linear part 
of the curve of growth. This is in agreement with the 
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findings in ICavrel et al.l (|2004ll , where the ctew is in¬ 
dependent of EW. A penalty is already incurred when 
trying to fit two lines when there is in fact only one. 


For our Gaussian line profiles, the equivalent width, 
EW, is given by EW = \P^aw, so that error propa¬ 
gation gives = 2ttw'^< 7 "^. This means that equa¬ 
tion da can also be interpreted as a relation between 
the uncertainties of the EW of a blended line to the 
corresponding value in isolation. 

We checked our result by Monte-Carlo (MC) simu¬ 
lation using two weak lines of the same strength with 
a = 6 = 0.1. Figure [5] illustrates the result of 1000 noise 
realizations. The MC test confirms that the ctew is in¬ 
dependent of EW (and therefore of line strength) for 
weak lines. For the MC run we chose a rather coarse 
(step size A A = FWHM/2) sampling, to show that the 
sampling is not particularly important for the valid¬ 
ity of the analytical result - as long as it is not overly 
coarse. We note that the way we set-up the MC run 
leads to correlated noise among the various S/N cases. 
This explains that shifts relative to the analytical result 
appear the same for all S/N. All events for a particular 
S/N are independent. To double gew in comparison to 
the case in isolation one needs a blending line as close 
as 


1.00 


s, 

¥ o-io 


0.01 


Fig. 6 Uncertainty gew of one component in a blend 
of two Gaussian components as a function of their sep¬ 
aration (both given in units of w). The analytical result 
according Eqs. dS]) is depicted by the dashed lines, re¬ 
sults of a Monte-Carlo simulation by solid dots. ±3 ct 
error bars (so small that they are barely visible) indi¬ 
cate the uncertainty of the Monte-Carlo result. Four 
cases with S/N=12.5, 25, 50, and 100 per pixel are 
shown from top to bottom. 


S/N=12.5 


1000 events, Sx=FWHM/2 


S/N=100 •' 


1 

Line separation [w] 


10 


4.3.2 Probing the line-to-line separation 
needed for an accurate chemical analysis of 
blended lines 


AAdoubie = \/2 ln(4/3) w « 0.759 • w 

K. 0.322-FWHM (7) 

At this separation two Gaussians appear essentially 
like a single one. If we instead look at the separation 
needed in order to derive abundances with an accu¬ 
racy of 0.1 dex or better, the necessary separation is 
^ 0.6 X FWHM. This value should be compared to the 
tests conducted and described in the next section. 

It might appear surprising that blended compo¬ 
nents can be separated to a distance of ~ 3/5 x FWHM 
according to the simple model presented in this section. 
However, the key property which makes this possible is 
the a-priori knowledge of the width of an isolated line. 
Close blends can be identified by a certain degree of 
additional broadening beyond what is expected for an 
isolated line. Conversely, the separation becomes more 
difficult when the intrinsic line width is not, or only 
approximately, known. 

In the next section we test this in practice, but we 
also change perspective and address the question of 
what happens if uncertainties in the wavelength cali¬ 
bration or measured radial velocity lead to changes of 
the position of a blend with respect to the expected po¬ 
sition based on the laboratory rest wavelengths of the 
blending components. 


Following the same assumptions as outlined in the pre¬ 
vious subsection we test the needed line-to-line separa¬ 
tion to recover the stellar abundances in practice. For 
this purpose we use synthetic spectra (mock data) that 
mimic those to be observed by 4MOST in the HRS 
setting. The two lines are weak, yet not so weak that 
they mistakenly could be taken as noise. The only free 
parameters are the two abundances, which we try to 
recover from spectra synthesized with different line-to- 
line separations (in units of FWHM) of the two blend¬ 
ing lines. The lines are treated in isolation and analysed 
without any ‘a priori’ (input) information from other 
lines of the same species. 

To test how close two blended lines can be, and 
remain useful for an abundance analysis, allowing for 
uncertainties of 0.1 dex, we start by selecting a region 
where this can be tested in the cleanest possible way. 

In the bluest part of the blue range, this will be hard 
to test, since many blends will arise simultaneously, es¬ 
pecially for the more metal-rich stars. Therefore, we 
choose a line that was clean in all the tests made for set¬ 
ting the wavelength region (see Sect. 14.2|) . The Cr I line 
at 427.48 nm provides a good test case, as it is clean in 
both dwarfs and giants from [Fe/H] = —3 to —1 dex. We 
then insert a line (here we chose Al as an example of a 
blend) in the core of the Cr line and in small steps move 
the centre of the Al line away from the centre of the Cr 
line. We generate synthetic spectra with separations in 
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steps of 0.1 XFWH1V0. The atomic data, adopted from 
a nearby Fe I line, was slightly altered to make this ‘Al’ 
line as stron^iil as the Cr line and visible at all three 
metallicities (—3,—2,-1) in both dwarfs and giants. To 
generalise this test to other elements that only show one 
blended line within the covered spectral range, we do 
not use abundance information from lines of the same 
species at other wavelengths. There are several other 
useful Cr and Al lines in the covered spectral range, 
but this is not the case for some of the other (heavier) 
elements. In this way we have two blending lines for 
which we have no a priori abundance information from 
other clean lines. 

The test was carried out using noise-free spectra 
generated with the spectrum synthesis code used in 
Sect. 14.21 where the lines were broadened to the afore¬ 
mentioned instrumental resolution. Only the abun¬ 
dances of the two blending lines were kept as free pa¬ 
rameters and the central wavelengths of these lines are 
assumed to be known. This allows us to simultaneously 
fit both lines and see if the abundances can be recov¬ 
ered with uncertainties better than ±0.1 dex. 



[Al/Fe] 


Fig. 7 Reduced for the automated simultane¬ 
ous fit to the isolated, blended lines of Al and Cr 
with a O.SxFWHM separation. Their abundances are 
shown on the x- and y-axis, respectively. Contours are 
smoothed with a kernel of 0.1 dex width, and the ’±’ 
marks the best solution (0.2,0.0) vs the true value ’dot’ 
(0.2,0.1). 

We set up a suite of 10 synthetic spectra with 
[Cr/Fe]= 0.1 and [Al/Fe]= 0.2 dex and with different 
line-to-line separation. These spectra were blindly anal- 

i.e., finer steps were used, compared to the step size in pre¬ 
vious section, to probe the separation needed. 

i.e., reach the same residual intensity 


ysed in three different ways by several co-authors. First 
we conduct an MC test on a randomly selected spec¬ 
trum of a metal-poor dwarf star (Teff/logg/[Fe/H]/Vt: 
5500K/4.0dex/-2.0dex/2kms-i). 

1) We create a grid of synthetic spectra with Al and 
Cr abundances varying in steps of 0.1 dex (the adopted 
allowed uncertainty) and derive the abundances using 

minimization. For this purpose we adopt a nominal 
uncertainty of 3% for the y^-fitting i.e., the uncertainty 
on the data to be fitted, which roughly corresponds 
to an overall 0.1 dex abundance uncertainty. At a sep¬ 
aration of 0.3xFWHM we recover the Al abundance 
([Al/Fe]= 0.2 dex) exactly, while we obtain a Cr abun¬ 
dance lower than the input, at [Cr/Fe]= 0.0 dex, which 
is at the limit of the allowed uncertainty. This case is 
shown in Fig. [7] for which we find a reduced of 0.93. 

2) A manual test was made on the same spec¬ 
tra by blindly analysing them in MOOG and treating 
them as real observed spectra for the same dwarf star 
(5500K/4.0dex/—2.0dex/2kms“^). At a separation of 
OxFWHM the input Cr and Al abundances could not 
be recovered, but we find several possible solutions. At 
a separation of O.SxFWHM or more, the input abun¬ 
dances can be recovered exactly. This agrees within 
±0.1FWHM with the value of 0.6 found in the previous 
theory section (14.3.11) . At 1 x FWHM the separation of 
the lines is visible, and the abundances can be derived 
as if the lines were not blended. 

In one test spectrum we introduced an extra lincP^ 
with an abundance of 0.1 dex into the blend as an un¬ 
known blend. This increased the derived Cr abundance 
by 0.1 dex. 

Finally, if we assume the abundance of one of the 
two blending lines to be known (e.g., Cr), based on 
measurements from other clean Cr lines, the abundance 
of the Al line can be perfectly recovered regardless of 
FWHM separation even when the centres of the two 
blending lines coincide. 

The observed 4MOST spectra will not necessarily 
be shifted perfectly to rest wavelength, but the radial 
velocity might be off by ^ lkms“^. According to an¬ 
other scientific requirement, the radial velocity must 
be accurate to within ±1 — 2kms“^. In metal-rich gi¬ 
ants most of the lines will be very strong or even satu¬ 
rated, which means that the line profiles will no longer 
maintain a Ganssian shapJ^. Gombining all these fac¬ 
tors and furthermore adding the radial velocity as a 
free parameter will aggravate the abundance recovery 
from a blend containing a saturated and a weaker line. 
Thus, by considering a short wavelength range consist¬ 
ing mainly of non-Gaussian shaped strong lines, could 
lead to problems where an uncertain radial velocity no 
longer causes a simple systematic offset. 

0.15 A from the Cr line 

This can also afffect the accuracy of the radial velocity shift, 
since saturated lines tend to yield more uncertain shifts. 


www.an-journal.org 


© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 































12 


C. J. Hansen et al.: Stellar science from a blue wavelength range 


3) Therefore a different test was carried out on a 
’worst case scenario’, namely a more metal-rich giant 
(4500K/2.0/—1.0/2kms“^) with numerous blends in 
the blue region. We used 0.1 dex lower abundance0 
for this test. The saturation plus offset radial velocity 
can hide the weaker line and return an abundance that 
is off by up to 0.3-0.4 dex (see Fig. [8l). To obtain the 
correct abundance a larger separation (> IxFWHM) 
is needed. The shift of the line centre (mock data off¬ 
set by 0, 0.5, and lkms“^see Fig. ED is barely visible 
in the broad line with an almost flat core of the satu¬ 
rated line. Only under careful (x^) inspection this will 
be detected. In an automated code this could result in 
a false minimum, that would overestimate the abun¬ 
dance of the weaker line by up to 0.3-0.4 dex, while the 
stronger line could be reproduced within ~ ±0.2 dex. 



Fig. 8 Mock data (observations - black diamonds) 
with different radial velocity shifts (i.e., offset by 0, 
0.5, 1 kms“^ from top to bottom) compared to synthe¬ 
ses with different [Cr/Fe] and [Al/Fe] abundances. All 
spectra have been broadened to the 4MOST resolution. 
Residuals (observations - syntheses) show the quality 
of the fit. At an offset of 0.5kms“^ the syntheses are 
equally good/bad. In order to derive abundances that 
are exact within ±0.1 dex the radial velocity shift needs 
to be exact (0kms“^) or the separation of the blending 
lines needs to be larger than IxFWHM. 


i.e., [Cr/Fe]=0.0 and [Al/Fe]=0.1 dex 


The last test (3) illustrates future limitations in the 
automated pipeline abundance analyses of ’worst case 
scenarios’ which is needed for a survey like 4MOST. 
Not surprisingly, the more metal-rich giants are the 
hardest cases from which to recover the abundances for 
each of the lines. The lines are extremely strong (sat¬ 
urated or even damped) in metal-rich giants, making 
the abundance recovery even harder. This is why the 
blue range with all its blends can generally not be used 
for more metal-rich disk and bulge stars. This does not 
mean that this range does not contain any useful lines. 
Section 14.21 shows that there are several clean lines in 
metal-rich giants for which we can derive abundances 
within the adopted uncertainty (±0.1 dex). Special at¬ 
tention needs to be paid to very metal-rich giants to 
ensure that the radial velocity shift is accurate since 
this could degrade the abundances derived from strong 
blended lines. 

Such an extreme case would be very unlikely at a 
metallicity lower than [Fe/H]< —1, where there are 
fewer (saturated) lines, and if the lines are not stud¬ 
ied in isolation, a ‘wrong’ radial velocity (±lkms“^) 
should reveal itself by offsetting the surrounding lines. 

This shows that as long as we deal with two weak, 
blended lines, these can be used for precise abun¬ 
dance measurements (±0.1 dex) in the 4MOST HRS 
spectra, provided the line centres are separated by 
^ 0.5xFWHM if all other parameters and input (e.g., 
atomic data, stellar parameters, and radial velocity) are 
well constrained. If we know the abundance for one of 
the blending elements from other clean lines, we can re¬ 
cover the abundance of the other component perfectly 
even if the two lines overlap. On the other hand, if the 
radial velocity is uncertain and the line is saturated or 
contains unknown blends, we will need a larger sepa¬ 
ration (> IxFWHM) to recover the abundances from 
the blending lines. 

5 Conclusions and outlook 

The future survey instrument 4MOST will provide data 
that are crucial for making advances in cosmology, ac¬ 
tive galactic nuclei, and the formation of galaxies. For 
the latter, the Milky Way will be used as a model or¬ 
ganism. 

4MOST will allow for studies on large scales to un¬ 
derstand galaxy formation as well as on small scales 
by tracing the origin of heavy elements. Both of these 
are among the open questions in astrophysics to date, 
and they cannot be answered to the same extent using 
green and red spectral studies only. In particular, the 
important science case of metal-poor stars cannot be 
studied and answered using only redder spectra. 

With the current design plan, the instrument will 
be able to simultaneously observe up to 2400 objects 
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in one pointing, using fibers feeding high- and low- 
resolution spectrographs. To be able to measure the 
abundances of at least the 20 elements that are most 
important for constraining the chemical evolution of 
the Galaxy, the blue arm of the 4MOST high-resolution 
spectrograph must provide a resolving power of i? ^ 
18 000-20 000, and a sampling of 2.5-3 pixels per reso¬ 
lution element. 

The blue spectral range is very crowded and there¬ 
fore many lines blend. Here we have shown that accord¬ 
ing to theory, the line centres of two weak, blending 
lines need to be separated by ~ O.OxFWHM to make 
an abundace recovery within ±0.1 dex uncertainty. Our 
subsequent tests on synthetic spectra showed that the 
separation needs to be ^ O.SxFWHM in order to re¬ 
cover accurate abundances (within ±0.1dex) from both 
components (unless we a priori know the abundance of 
one of the elements). The precision of the radial veloc¬ 
ity determination plays an important role in how accu¬ 
rately we can recover the abundances from two blended 
strong lines. 

To date only three high-resolution instruments are 
capable of observing all the way down to the atmo¬ 
spheric cut-off while most others stop ^ 50 — 100 nm 
redward of this. However, the 300-400 nm region is im¬ 
portant when observing halo stars no matter if the sci¬ 
ence cases focus on constraining nuclear formation pro¬ 
cesses of elements with Z> 38 or searching for the ori¬ 
gin of carbon enhanced metal-poor stars, which could 
be some of the first stars formed after the Population 
HI stars. However, to unravel the origin of stars, nu¬ 
clear processes, or large scale formation scenarios of 
the Milky Way, a large sample of stars is needed at 
[Fe/H] < —3. The efficiency of the fibres needs to im¬ 
prove to cover a wider wavelength range at higher res¬ 
olution in the blue to make best use of these blue-most 
wavelength regions. Alternatively, a stronger initiative 
in the community focusing on blue-to-near-UV science 
is needed, and even more time needs to be allocated 
to high-resolution slit spectroscopy covering a deeper 
follow-up of these blue wavelength ranges. 
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